Determination of Inter-Phase Line Tension in Langmuir Films 
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A Langmuir film is a molecularly thin film on the surface of a fluid; we study the evolution of a 
Langmuir film with two co-existing fluid phases driven by an inter-phase line tension and damped 
by the viscous drag of the underlying subfluid. Experimentally, we study an 8CB Langmuir film 
via digitally-imaged Brewster Angle Microscopy (BAM) in a four-roll mill setup which applies a 
transient strain and images the response. When a compact domain is stretched by the imposed 
strain, it first assumes a bola shape with two tear-drop shaped reservoirs connected by a thin tether 
which then slowly relaxes to a circular domain which minimizes the interfacial energy of the system. 
We process the digital images of the experiment to extract the domain shapes. We then use one of 
these shapes as an initial condition for the numerical solution of a boundary-integral model of the 
underlying hydrodynamics and compare the subsequent images of the experiment to the numerical 
simulation. The numerical evolutions first verify that our hydrodynamical model can reproduce the 
observed dynamics. They also allow us to deduce the magnitude of the line tension in the system, 
often to within 1%. We find line tensions in the range of 200-600 pN; we hypothesize that this 
variation is due to differences in the layer depths of the 8CB fluid phases. 

PACS numbers: 68.18.-g, 68.03.Cd, 61.30.Hn 



Line tension, the two dimensional analog of surface 
tension, is the free energy per unit length associated 
with the boundary between two phases on a surface. 
In this paper we explore a method for measuring the 
inter-phase line tension in Langmuir layers, the quasi- 
two-dimensional surface layers of polymers, lipids or liq- 
uid crystals that exist at gas-liquid and liquid-liquid in- 
terfaces. Langmuir layers often separate into multiple 
domains signaling the coexistence of different phases [l| . 
The boundaries of such domains are curved, yielding a 
line force per unit length normal to the phase boundary 
and tangent to the surface containing the Langmuir layer 
with a magnitude that is the product of the line tension 
and the curvature of the interphase boundary. 

Attempts to measure the line tension in various sys- 
tems have multiplied over recent years. One motivation 
is to better understand the forces which govern the shape 
and influence the function of biological membranes; cell 
membranes consist of a mixture of cholesterol, lipids, and 
proteins that can form domains with various structures 
and functions. Model membranes, including supported 
bilayers Q , vesicles Q and Langmuir monolayers [l| show 
macroscopic phase separation, with geometry driven by 
line tension. 

Line tension between fluid Langmuir phases has most 
often been measured by watching the relaxation of 
stretched domains toward an energy-minimizing circu- 
lar shape. The relaxation of large perturbations, such 
as bola-shaped domains (two teardrop-shaped resevoirs 
tethered together with a line of nearly constant thickness) 



have been modelled only heuristically; models to extract 
line tension [1, [f| [|| approximated the bola shape as two 
perfectly round discs connected by an infinitesimally thin 
tether, which is far from the true form. The dynamics 
of linearized perturbations of circular domains are better 
understood [!, 0, but these perturbations are diffi- 
cult to measure accurately in the small amplitude limit 
where they obtain validity. Due to these problems, the 
error bounds of previous line tension measurements have 
been no better than ±20%. 



Our group recently developed a manageable model Q 
of the experimentally observed relaxation dynamics of 
two fluid phases within a Langmuir film. The model is 
both analytically tractable and allows an efficient, ac- 
curate and stable numerical solution via a boundary- 
integral technique. 



In this article we directly compare the numerical re- 
sults of our model to experimental results on a Lang- 
muir layer with two fluid phases corresponding to differ- 
ent multilayer thicknesses, and test both the validity of 
the model and the precision of the line tension measure- 
ments resulting from that comparison. We expect this to 
set the stage for further accurate and precise studies of 
line tension as a function of temperature, composition, 
and other variables. 



I. EXPERIMENTAL 

We conduct our measurements on Langmuir films com- 
prised of 4'-8-alkyl[f ,f'-biphenyl]-4-carbonitrile (8GB) 
deposited on a subfluid of pure water. The 8CB exists as 
a smectic liquid crystal with stacked molecular bilayers 
on top a simple monolayer at the water surface 0, 0, fiol ] . 
Consequently, multiple phases each consisting of a differ- 
ent odd number of layers (i.e. monolayers, trilayers, etc.) 
can simultaneously exist within the film. 

Relaxation in Langmuir layers is driven by intermolec- 
ular forces between the surface molecules and also be- 
tween the layer and the subfluid. In some systems, elec- 
trostatic forces in the Langmuir layer (primarily dipole- 
dipole repulsion) drive interesting pattern formation such 
as circle-to-dogbone transitions and labyrinth formation 
[Hi 13- We choose the 8GB multilayer system consid- 
ered here because the electrostatic effects are probably 
negligible, in that a symmetric bilayer is added at each 
step. No jump in surface potential, which determines 
the effective dipole moment density, is observed after the 
triple layer. 

Consequently, in this system the intermolecular forces 
are well-modeled as a line tension at phase boundaries, 
which causes the film to coalesce into spatially-distinct 
phase-domains. Any domain strained into a non-circular 
shape will relax to the energy- minimizing circular config- 
uration, driven by the inter-phase line tension. 

The 8CB forms a smectic phase at the water surface 
which behaves like a two-dimensional fluid. The surface 
viscosity can be estimated from the bulk viscosities of 
the smec tic p hase; it is less than 100 times the viscosity 
of water jl3l . [l4| , so that for domains we consider with 
thickness less than 100 nm, the surface viscosity is negli- 
gible as long as the domain size is ^> 10/im. 

The 8CB (Sigma- Aldrich, 98% pure) is further purified 
by chromatography. We dissolve 8GB in hexane (Fischer, 
Optima grade) spreading solution, which is deposited on 
the surface of water (PureLab+ system, and passes the 
shake test) in a clean trough (mini-trough, KSV). After 
deposition, the hexane evaporates, leaving an 8CB layer 
on the water surface. The trough has a pair of movable 
barriers to change the water surface area available to the 
8GB film and thereby control the surface pressure. At 
room temperature, surface pressures ^6.5 mN/m pro- 
duce a stable coexistence of a tri-layer over the entire 
surface interspersed with thicker domains [§, [To| . We 
image the Langmuir layer using a homemade Brewster 
Angle Microscope (BAM) [IH, [l6|, [l?]], which produces 
grey-scale images showing more thickly-stacked domains 
in brighter shades against the dark, thinly-stacked back- 
ground. 

We stretch the domains by shearing the subfluid and 
then use the BAM to observe the subsequent relaxation, 
which is recorded on a computer at 30 frames per second. 
To shear the subfluid we use a 4-roll mill [l8|, [l|| H(| con- 
trolled by a stepper motor. The rolls are made of black 
Dehin which is hydrophilic and has no measurable effect 



2 




FIG. 1: Schematic of the four-roll-mill flow profile. The ge- 
ometry of the rollers, a — 6.6 mm and b = 10.5 mm, produces 
maximum homogeneity in the extension rate (l8lll9. 20],. The 
shaded ellipse shows the area illuminated by the laser beam. 

on surface pressure. We adjust the water level to be ex- 
actly the same height as the upper edges of the rolls in 
order to minimize the distortion of the fluid surface re- 
sulting from contact with the rolls. As shown in Figure 
[TJ the 4-roll mill provides symmetric shear forces about a 
central stagnation point on the surface. This allows us to 
stretch a domain located at the stagnation point without 
imparting a net velocity to the domain and moving it out 
of the BAM's field of view. A controlled air stream ma- 
neuvers a domain into proper location at the stagnation 
point. Once the domain is in position, we activate the 
4-roll mill, and the domain stretches out, assuming the 
characteristic bola shape . Generally, we run the 4-roll 
mill at speeds of ~0.2 revolutions per second for about 5 
seconds. In our experiments, the Reynolds number plj 
of the flow during shearing is ~16. Because of the iner- 
tia in the subfluid, the domain continues stretching for 
several seconds after the mill has been stopped. 

II. HYDRODYNAMICS 

Our model Q describes the dynamics of a Langmuir 
layer consisting of two phases: an isolated phase-domain, 
f2, of finite area surrounded by a second Langmuir phase, 
il c , which extends infinitely in the horizontal direction. 
The Langmuir layer is modeled as a flat, two-dimensional 
fluid. We assume that the subfluid is infinitely deep. 
Both the Langmuir layer and the subfluid are assumed 
to be incompressible on the timescale of the relaxation 
experiments. Thus, the Langmuir domain Q will have a 
fixed area, A*. 

For the Langmuir layer, the incompressibility condition 
in relaxation driven by line tension alone corresponds to 
a Gibbs elasticity Kn 3> A/VA*, where A is the inter- 
phase line tension [22j . For the experimentally accessible 
range we estimate that V 'A* > 10/im and conservatively 
choose an upper limit on the line tension of A < InN, 
which yields an upper bound on the Gibbs elasticity, 
Kg ^S> 10~ 4 iV, for significant compressibility. Thus, al- 
most any Langmuir layer liquid and many gasses will act 
as incompressible fluids in line-tension driven relaxation 
processes (see [22J and references therein). In the spe- 
cial case of multilayers considered here, it is additionally 
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conceivable that the number of layers might change dur- 
ing the relaxation process, leading to an area change. In 
practice, the thickness of the layers and the area of the 
domains were observed to be constant for the duration 
of the experiments reported here. 

As we discuss above, dimensional analysis indicates 
that for 8CB the energy dissipated by viscous shearing 
within the Langmuir layer is much less than the amount 
dissipated by viscous shearing of the subfluid; we there- 
fore model the Langmuir layer as inviscid. Furthermore, 
the subfluid can be treated in the Stokesian limit, where 
its inertia is negligible. 

We non-dimensionalize the dynamics in terms of a 
characteristic length, time, and mass 
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r/A* 
A 



M* = r]'L*T* 



(1) 



respectively; here A is the inter-phase line tension and rf 
is the subfluid viscosity. Essentially, the relaxation of the 
domain f2 is driven forward by the line tension between 
phases, and slowed by the viscosity of the subfluid. 

The model ultimately yields an equation of motion for 
the boundary curve, T, separating f2 and il c . As the 
boundary is isotropic, it suffices to determine the nor- 
mal velocity, U, to specify the domain's evolution. We 
obtain i| 



U 



ds 



(2) 



where fi is the outward unit normal vector to T, s is the 
arclength measured in a right-handed sense, and ^>(s) is 
the velocity streamfunction restricted to the boundary of 
the domain. This is computed as a boundary integral, 
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k(s') i(s') ■ Q(s, s) 



ds' 



(3) 



where k is the curvature, t is the unit tangent vector, and 
Q(s,s') is a unit vector pointing from T(s) to r(s'). A 
derivation and discussion of this formulation is given in 

We implement a numerical solution in MATLAB. The 
problem is extremely stiff numerically; explicit integra- 
tion methods are very susceptible to high-wavenumber 
instabilities. This can be ameliorated by operator split- 
ting, following the ideas of Hou et al. [23j |, While such 
a splitting in not immediately apparent in the formula- 
tion above, the formulation in Lubensky and Goldstein 
[24j and Heinig et al. [25| can be used to show that the 
high-frequency modes of T are asymptotically governed 
by a much simpler evolution law, namely motion by mean 
curvature. 

As in [HI , using an intrinsic description of the bound- 
ary allows an accurate implicit solution for the high- 
wavenumber modes, avoiding numerical instabilities. We 
represent the boundary with an equal-arclength dis- 
cretization. Derivatives are computed pseudo-spectrally 



(261 |27I | , and the boundary integral is computed using the 
16-panel closed Newton-Cotes formula which guarantees 
high-order spatial accuracy. It is straightforward to solve 
the evolution by mean curvature implicitly and to high 
accuracy [23j |. We proceed by using Strang splitting [28| 
with the mean-curvature step implemented implicitly and 
the boundary integral velocity minus the mean curvature 
velocity computed explicitly. 

Numerically, we observe a slow drift of the grid which 
forces us to regularly correct the arclength discretization 
— this is done using spectral interpolation and Newton- 
Raphson iteration. Also, it is necessary to filter the 
highest-frequency modes of T (whose numerical accuracy 
is poor due to the discretization anyway); we convolve 
the spectrum with a smooth filter and retain roughly 
two-thirds the spectrum. Details of the numerical imple- 
mentation are available in 1 29] . 



III. FINDING THE BOUNDARY CURVES 

To analyze a set of experimental photographs we must 
first determine the location of the boundary curve in each 
one. A grey-scale photograph is a map from each pixel 
(x, y) to the brightness of the image at that location, 
B(x,y). The edge of the domain is located in the region 
of rapid transition from black to white, where ||V-B|| is 
large. We compute V.B and ||VB||, using code devel- 
oped by Fisher et al [3(|. We execute a curve-tracing 
algorithm that "walks" around the edge of the domain, 
staying in the thin region where ||V-B|| is large. As the 
algorithm traverses the boundary it marks points, which 
we subsequently use as a discrete representation of the 
boundary curve. 

The placement of the edge can be quickly verified vi- 
sually. We also have a quantitative check at our disposal. 
The domain area is conserved; if the edge is placed too 
far to the outside or inside then as the perimeter of the 
domain decreases during relaxation the computed area 
of the domain decreases or increases, respectively. This 
relationship allows us to calculate the (average) distance 
by which the edge is displaced in the normal direction. 
In most data sets we see no correlation between perime- 
ter length and domain area. When this effect is seen, the 
implied displacement is never more than two pixels, and 
we can move the edge in the normal direction to correct 
for the displacement (this was done in data sets B and E 
reported below). 

The greatest obstacle to determining the precise loca- 
tion of the edge is diffraction, which blurs the edge and 
produces a bright ring of constructive interference. Al- 
though the average normal displacement of the curve is 
very small, diffraction may cause the edge to be off by 
several pixels locally (this problem affects both human 
visual perception and computer algorithms). 
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FIG. 2: A series of snapshots (for data set B) showing 
the relaxation of a stretched domain to a circle. The snap- 
shots are separated by 2.85 seconds w 0.558 T*. The first 
photo is superimposed with a curve marking the bound- 
ary determined by our edge-finding algorithm; this is the 
initial condition for the simulation of the domain relax- 
ation. In subsequent photos the superimposed curve shows 
the simulated shape of the bola, based on that initial con- 
dition. The simulation accurately reproduces the observed 
dynamics even for complex, asymmetrical domains. The 
fine tension of this domain is found to be 390 ± 3 pN, 
where the uncertainty corresponds to the standard devia- 
tion of T*. A movie of this relaxation is available from 
http : //www . math . hmc . edu/~ajb/bola/ simulation.mpeg. 



IV. DETERMINING A VIA 
EXPERIMENT/SIMULATION COMPARISON 

Our equation of motion (|2|) is written in time units of 

the characteristic relaxation time, defined by T* = 71 ^* . 
We calculate the line tension, A, given these other values. 
The domain's area, A*, is determined from the photo- 
graph once the boundary curve is found, and the subfiuid 
viscosity, rj' , is estimated from tabulated values (adjust- 
ing for temperature). We determine T* by simulating the 
evolution and matching timescales between the observed 
relaxation (snapshot times are recorded in seconds) and 
the simulated relaxation (done in units of 2 1 *) . 

We choose an initial condition taken from one of the 
snapshots and simulate the subsequent relaxation. For 
each photograph after this first one, we match timescales 
by finding the time in the simulation when the shape of 
the simulated domain most closely matches the domain 
shape in the photo. We search through the discrete time- 
steps of the simulation and compare the shape at each 
step to the shape taken from the photo. Each snapshot 
gives us a value for T*, computed as 

where tj and to are the observed times of the compari- 
son snapshot and initial-condition snapshot, respectively, 
and Tbest is the time which elapses in the simulation be- 
tween T = (the initial condition) and the time at which 
the best-matching shape occurs. 

To measure how closely two domains (i.e. image- 
processed experiment and numerical simulation) match 
we use the Symmetric Difference Metric (SDM), which 
is determined by overlaying the domains and comput- 
ing the total area which lies in one or the other but not 
both. For each photo, we search through the time-steps 
of the simulation to find the step at which the SDM be- 
tween simulation domain and the experimental domain 
is minimized; Figure [3] provides a visual illustration of 
this process. The minimum SDM over the simulation 
provides a measure of how well each photo matches some 
shape which occurs in the simulation] we expect the same 
value for T» from every photo. We compute the mean 
value of the set of T* 's from all the photos to determine 
the line tension, and the standard deviation of these TVs 
provides an estimate for the precision in the resulting 
measurement of A. 

To simulate the relaxation we must know the compo- 
nent of the subfiuid velocity which exists independent of 
(i.e. is not directly produced by) the relaxation. Un- 
fortunately, we cannot directly measure the subfiuid ve- 
locity. Instead, we choose an initial snapshot when the 
subfiuid is relatively quiescent and run the simulation un- 
der the approximation that the "independent" subfiuid 
velocity is zero. We find, however, that the violation 
of this approximation is one of the largest sources of sys- 
tematic error contributing to mismatch between observed 
and simulated domain shapes. We could choose a later 
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FIG. 3: Three images showing the same photo overlayed 
with simulated boundary curves from three different times in 
the simulation. The center image shows the simulated curve 
at the time, T best , when it most closely matches the photo; 
the left and right images show the simulated curve at T — 
T best - 0.1T, and T = T best + 0.1T», respectively. 



initial condition, waiting until remnant subfluid velocity 
is negligible; however, this means throwing out a large 
portion of the data (often all of it). 

The type of motion which persists in the subfluid for 
the longest time is solid-body motion; other types of 
motion are viscously damped. We therefore correct for 
solid-body motion in the post-simulation timescale fit- 
ting. Whenever we compare two shapes, we do not di- 
rectly compute the SDM between them, but instead de- 
termine the minimum SDM which can be achieved by 
positioning one on top of the other using a solid-body 
motion. This greatly reduces the SDM and allows us to 
achieve excellent fits for data which would otherwise be 
rendered worthless by remnant subfluid velocity. 

Following Q , we also measure the line tension by mea- 
suring the relaxation of small elliptical deformations of 
the boundary in the near circular limit, which we refer 
to as X e . The snapshots of the domain boundary are im- 
age processed, and FFT techniques are used to extract 
the amplitude of the elliptical (n = 2) deformation. We 
then fit the exponential relaxation rate of this mode in 
the small amplitude limit. 

Comparison data for six separate relaxations is pre- 
sented in Table HI In six time series of different domains, 
the mean SDM between the experimental and simulated 
domains was 1.5-4% of the domain's area, indicating that 
the proposed hydrodynamical model of domain evolution 
reproduces the shapes quite well; this is clear from the 
comparison snapshots in the evolution in Figure The 
areas of the domains were constant across the time series 
to 0.1-1%, well within the uncertainties of the measure- 
ment due to diffraction at the domain edges. By match- 
ing time scales between the experiment and simulation, 
each photo after the first yields a value for in seconds; 
we deduce X avg from these values. The percentage de- 
viation of the values for T* from a set of photos ranges 
from 0.4%-4%, which also provide an error estimate on 
the line tension. 

The greatest variances in T» and the largest SDM val- 



Data Set 


Xavg [pN] 


0"A/ Xavg 


Average SDM 


0~Al A aV g 


Xe [ P N] 


A 


538 


3.4% 


3.9% 


1.0% 


468 


B 


390 


0.8% 


2.9% 


0.7% 


375 


C 


357 


1.7% 


3.0% 


0.3% 


362 


D 


570 


2.0% 


1.5% 


0.1% 


606 


E 


479 


4.0% 


1.5% 


1.0% 


485 


F 


191 


0.4% 


3.4% 


0.5% 


217 



TABLE I: Line tension values and error estimates for six data 
sets. Here Xavg is the line tension values averaged over snap- 
shots from large aspect ratio of bolas through the relaxation 
to nearly circular domains. The percentage error in this mea- 
surement is estimated by normalizing the standard deviation 
of these measurements by the average line tension, a\/X a vg- 
The average of the SDM (Symmetric Difference Metric) nor- 
malized by the domain area in these snapshots reflects how 
closely the numerical simulation reproduces the experimental 
results. The standard deviation of the domain area, a a di- 
vided by the average area, A avg , is a proxy for the error in 
the image processing of the boundary. Finally, A e is the line 
tension estimated from small amplitude perturbations from 
the final circular domain shape. 

ues (i.e. shape mismatches) occurred in those data sets 
where either (1) other domains nearly touched the do- 
main of interest or (2) remnant subfluid shearing was 
particularly problematic. Provided that reasonably iso- 
lated domains can be produced and the subfluid flow can 
be well-controlled, it is possible to determine the line ten- 
sion to a precision of ~1%. 

Finally, we note that the line tension estimates, A e , 
from small perturbations of the final circular shape are 
off by up to 13 % - this is comforting, in that it is consis- 
tent with our results and variations observed in previous 
work Q. It suggests that our new methodology is both 
accurate and precise. 



V. CONCLUSIONS 

In this paper we have described a method for determin- 
ing the line tension driving the evolution of Langmuir lay- 
ers. We are able to verify that our hydrodynamic model 
is consistent with the experiments and to determine the 
line tension with errors as small as 1%, more than an 
order of magnitude better than previous efforts. 

While we believe our measurement are accurate, it is 
striking that we have observed a wide variation in line 
tensions (191-570 pN) for the 8CB system. One fac- 
tor that contributes significantly to this variation is the 
thickness of both the compact domain and its surround- 
ings. Experiments reported elsewhere [32^ . using the re- 
laxation of small deformations generated with a different 
deformation technique, systematically explore the depen- 
dence on the thickness of a compact domain 1-15 bilay- 
ers thick (on top of an unpaired monolayer) in a trilaycr 
background. These experiments lead to a line tension, 
reflecting the elastic energy of the dislocation at the do- 
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main boundary, proportional to the Burgers vector [32|. 
In the experiments reported here, we estimate, from the 
observed brightness of the domains, that the lighter com- 
pact domains range from 10-24 bilayers thick and that the 
darker surrounding region is either three or five layers 
thick. Because of light scattering from the four-roll mill 
in our study the thickness of the dark region, expected to 
contribute to the line tension [3lj . is particularly uncer- 
tain. Other possible factors influencing the line tension 
include contamination and splitting of the boundary into 
two dislocation lines for very thick films [IJ. Our ex- 
pectation is that by quantifying this system and others 
more carefully we will be able to determine the dominant 
causes of line tension and generate reproducible results. 
This is a promising area for future investigation. 

The present model assumes that bulk viscosity dom- 
inates the relaxation and that both slip between layers 
and electrostatic effects are negligible. These conditions 
must be evaluated on a case-to-case basis. However, the 
current method can determine the line tension with any 
technique exploiting domain hydrodynamic response, in- 
cluding relaxation after coalescence of two domains [33[ 
or after stretching a domain with lasers tweezers [34j . It 
can also be generalized to more complex situations in- 
volving three-phase contacts. 

The boundary integral formulation here can be ex- 



tended to incorporate more general potential forces such 
as electrostatics (cf . [U [2f| ) and as such we believe that 
we have developed a valuable tool for deducing and veri- 
fying the form of the intermolecular potential in systems 
that exhibit more complex morphology such as circle- 
to-dogbone transitions and labyrinth formation [ill, E3] • 
This promises to be fertile ground for future research. 
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